Method and apparatus for I/Q mismatch calibration of transmitter

ABSTRACT

A method for I/Q mismatch calibration of a transmitter. First, a discrete-time signal is generated. Next, a corrected signal based on the discrete-time signal and a set of correction parameters A p  and B p  is obtained. Next, the corrected signal is converted to an analog corrected signal. Next, I/Q modulation is applied to the analog corrected signal and outputs a modulated signal. Next, a first, second, and third desired component measures and a first, second, and third image component measures with the first, second and third sets of the correction parameters A p  and B p  are respectively obtained from the modulated signal. Next, a fourth and fifth set of correction parameters A p  and B p  are obtained based on the first, the second, and the third desired component measures as well as the first, the second, and the third image component measures. Next, a fourth desired component measure and a fourth image component measure with the fourth set of correction parameters A p  and B p  and a fifth desired component measure and a fifth image component measure with the fifth set of correction parameters A p  and B p  are obtained from the modulated signal. Finally, a final set of the correction parameters A p  and B p  are selected from the fourth and fifth sets of correction parameters.

BACKGROUND OF THE INVENTION

The present invention relates to compensation of physical layer impairments on transmitter systems and particularly to a method and apparatus for I/Q mismatch calibration in transmitter systems.

DESCRIPTION OF THE RELATED ART

OFDM is a multi-channel modulation system employing Frequency Division Multiplexing (FDM) of orthogonal sub-carriers, each modulating a low bit-rate digital stream. The simplest way to describe an orthogonal frequency-division multiplexing (OFDM) signal is as a set of closely spaced frequency-division multiplexed carriers. While this is a good starting point for those unfamiliar with the technology, it falls short as a model for analyzing the effects of signal impairment.

The reason it falls short is that the carriers are more than closely spaced; they are heavily overlapped. In a perfect OFDM signal, the orthogonality property prevents interference between overlapping carriers. This is different from the FDM systems. In FDM systems, any overlap in the spectrums of adjacent signals will result in interference. In OFDM systems, the carriers will interfere with each other only if there is a loss of orthogonality. So long as orthogonality can be maintained, the carriers can be heavily overlapped, allowing increased spectral efficiency.

Table 1 lists a variety of common analog signal impairments and their effects on both OFDM signals and the more familiar single-carrier modulations such as quadrature phase-shift keying (QPSK) or 64-QAM (quadrature amplitude modulation). Most of these impairments can occur in either the transmitter or the receiver.

Impairment OFDM QPSK I/Q gain balance State spreading Distortion of (uniform/carrier) constellation I/Q quadrature skew State spreading Distortion of (uniform/carrier) constellation I/Q channel mismatch State spreading State spreading (non-uniform/carrier) Uncompensated State spreading Spinning constellation frequency error Phase noise State spreading Constellation phase (uniform/carrier) arcing Nonlinear distortion State spreading State spreading Linear distortion Usually no effect State spreading if not (equalized) equalized Carrier leakage Offset constellation Offset constellation for center carrier only (if used) Frequency error State spreading Constellation phase arcing Amplifier droop Radial constellation Radial constellation distortion distortion Spurious State spreading or State spreading, shifting of affected generally circular sub-carrier

For cost reasons, analog in-phase and quadrature (I/Q) modulators and demodulators are often used in transceivers especially for wide bandwidth signals. Being analog, these I/Q modulators and demodulators usually have imperfections that result in an imperfect match between the two baseband analog signals, I and Q, which represent the complex carrier. For example, gain mismatch might cause the I signal to be slightly different from the Q. In a single-carrier modulation system, this results in a visible distortion in the constellation—the square constellation of a 64-QAM signal would become rectangular.

To better understand how gain imbalance will affect an OFDM signal, look at the equations describing each individual sub-carrier. In the following analysis, it's important to keep in mind that, while an individual sub-carrier is analyzed, the I/Q gain imbalance error is on the signal that is the composite of all sub-carriers.

In the equation (1), C_(k,m) is a complex number representing the location of the symbol within the constellation for the kth sub-carrier at the mth symbol time. For example, if sub-carrier k is binary-phase-shift-keying (BPSK) modulated, then C_(k,m) might take on values of ±1+j0. The complex exponential portion of equation (1) represents the kth sub-carrier, which is amplitude- and phase-modulated by the symbol C_(k,m). Therefore: C_(k,m)(e^(j2πkΔft))  (1) where Δf is the frequency space between the sub-carrier. Using Euler's relation, the equation (1) can be rewritten as: C_(k,m)(cos(2π·kΔft)+j sin(2π·kΔft))  (2)

Now add a term “β” to represent gain imbalance. For a perfect signal, set β=0. As shown, the gain imbalance term will also produce a gain change. This was done to simplify the analysis. Therefore: C_(k,m)((1+β)cos(2π·kΔft)+j sin(2π·kΔft)  (3)

The equation can be rearranged and this can be rewritten as the sum of a perfect signal and an error signal: C_(k,m)(cos(2π·kΔft)+j sin(2π·kΔft))+C_(k,m)β cos(2π·kΔft)  (4)

Finally, converting back into complex exponential notation, we get:

$\begin{matrix} {{C_{k,m}{\mathbb{e}}^{j\; 2\;\pi\; k\;\Delta\;{ft}}} + {\left( {C_{k,m}\frac{\beta}{2}} \right) \cdot \left( {{\mathbb{e}}^{j\; 2\;\pi\; k\;\Delta\;{ft}} + {\mathbb{e}}^{{- j}\; 2\;\pi\; k\;\Delta\;{ft}}} \right)}} & (5) \end{matrix}$

In words, the equation (5) shows that a gain imbalance produces two error terms. The first error term is at the frequency of the kth sub-carrier. The second error term is at the frequency of the −kth sub-carrier. The phase and magnitude of the error terms are proportional to the symbol being transmitted on the kth sub-carrier. Another way of saying this is that I/Q gain imbalance will result in each sub-carrier being interfered with by its frequency mirror-image sub-carrier. Persons skilled in the art will instantly recognize this as imperfect sideband cancellation.

The equation (5) has several implications. First, it is generally true that for sub-carriers used to carry data (as opposed to pilots), the symbol being transmitted at any given time on the kth sub-carrier is uncorrelated to the symbol on the −kth sub-carrier.

For a given sub-carrier, the lack of correlation from the mirror-image sub-carrier implies a certain randomness to the error. This results in a spreading of the sub-carrier's constellation states in a noise-like fashion. This is especially true for higher-order modulations such as 64-QAM. For lower-order modulations, such as BPSK, the error term from the mirror-image carrier has fewer states.

This can result in constellations where the BPSK pilot carriers of an 802.11a signal exhibit spreading that does not appear noise-like. Also, as the BPSK pilots do not have an imaginary component; the error terms associated with the pilot sub-carriers are real—so the spreading is only along the real (I) axis. Note that the phase relationships between the pilot carriers in an 802.11a system are highly correlated, so the errors introduced by quadrature errors are not random.

Quadrature skew produces error terms similar to those produced by gain imbalance. Quadrature skew occurs when the two oscillators used in an I/Q modulator or demodulator do not differ by exactly 90°. For a small angular error, it can be shown that the resulting error is orthogonal to the data. This is indicated by the j in front of the error terms in the equation (6). As with gain imbalance, the error generates energy at the kth and −kth sub-carriers. Again, the 802.11a BPSK pilots do not have an imaginary component, so the error term, which is now orthogonal, causes spreading along the Q axis. For the QPSK carriers in this example, the error is also orthogonal. However, unlike BPSK, a QPSK constellation doesn't look any different when rotated by 90°. (See the equation (6)):

$\begin{matrix} {{C_{k,m}{\mathbb{e}}^{j\; 2\;\pi\; k\;\Delta\;{ft}}} + {j{\frac{C_{k,m}\phi}{2} \cdot \left( {{\mathbb{e}}^{j\; 2\;\pi\; k\;\Delta\;{ft}} + {\mathbb{e}}^{{- j}\; 2\;\pi\; k\;\Delta\;{ft}}} \right)}}} & (6) \end{matrix}$

In both 802.11a and Hiperlan2, a channel estimation sequence is transmitted at the beginning of a burst. This special sequence is used to train the receiver's equalizer. The intended function of the equalizer is to compensate the received signal for multi-path distortion a linear impairment in the signal that is the result of multiple signal paths between the transmitter and the receiver. As the ideal channel estimation sequence is known by the receiver, the receiver can observe the effects of the channel on the transmitted signal and compute a set of equalizer coefficients.

In the transmitter, the channel estimation sequence is created by BPSK modulating all 52 carriers for a portion of the preamble. Not coincidentally, the equalizer consists of 52 complex coefficients—one for each sub-carrier. It should come as no surprise that each sub-carrier in the channel estimation sequence has the greatest influence on the equalizer coefficient computed for that same sub-carrier.

The channel estimation sequence, and the receiver algorithms that compute the equalizer coefficients, are not immune from signal impairments. Consider, for example, the effect of I/Q gain imbalance on sub-carriers +26 and −26 of the channel estimation sequence. Recall from equation (5) that each sub-carrier has two error terms: one at the same frequency as the sub-carrier, and one at the mirror image frequency. The I/Q gain imbalance will cause mutual interference between sub-carriers +26 and −26.

From the IEEE 802.11a standard, the sub-carrier modulation for the channel estimation sequence is defined to be C⁻²⁶=1+j0 and C₊₂₆=1+j0. Using these values in equation (5), one can easily determine that the two sub-carriers, when combined with the resulting error terms, will suffer an increase in amplitude. The equalizer algorithm will be unable to differentiate the error from the actual channel response, and will interpret this as a channel with too much gain at these two sub-carrier frequencies. The equalizer will incorrectly attempt to compensate by reducing the gain on these sub-carriers for subsequent data symbols.

The result will be different for other sub-carrier pairs, depending on the BPSK channel estimation symbols assigned to each.

With QPSK sub-carriers, the equalizer error caused by gain imbalance, or quadrature skew, results in seven groupings in each corner. Each QPSK sub-carrier suffers from QPSK interference from its mirror image. This results in a spreading to four constellation points in each corner. Each QPSK sub-carrier also suffers from a bi-level gain error introduced by the equalizer. This would produce eight groupings, except that the gain error is such that corners of the groupings overlap at the ideal corner state. Only seven groupings are visible.

I/Q Channel Mismatch

When the frequency response of the baseband I and Q channel signal paths are different, an I/Q channel mismatch exists. I/Q channel mismatch can be modeled as a sub-carrier-dependent gain imbalance and quadrature skew. I/Q gain imbalance and quadrature skew, as described above, are simply a degenerate form of I/Q channel mismatch in which the mismatch is constant over all sub-carriers. Think of channel mismatch as gain imbalance and quadrature skew as a function of a sub-carrier. It is still generally true that channel mismatch causes interaction between the kth and −kth sub-carriers, but that the magnitude of the impairment could differ between the kth and the (k+n)th carriers.

In order to eliminate the effects of the previously described impairments on the OFDM systems, various kinds of compensation circuits and methods have been proposed.

U.S. Application Publication No. 20020015450 discloses a method and an arrangement for determining correction parameters used for correcting phase and amplitude imbalance of an I/Q modulator in a transmitter. The transmitter includes an I/Q modulator and a corrector for correcting the phase and amplitude imbalance caused by the I/Q modulator. The arrangement has means for sampling the I/Q-modulated test signal to be transmitted, means for A/D-converting the signal samples taken from the test signal, means for I/Q-demodulating the signal samples digitally into I- and Q-feedback signals, means for determining the phase and amplitude imbalance caused by the I/Q modulator on the basis of the I- and Q-feedback signals, and means for determining the correction parameters of phase and amplitude on the basis of the determined phase and amplitude imbalance.

U.S. Pat. No. 6,208,698 discloses a quadrature modulator controlled by an imbalance estimator, as shown in FIG. 1. The estimator is entirely analog and includes two mixers 120, 121 to which are applied two carrier signals P1 and P2 derived from a carrier signal P obtained from a local oscillator and from a 90° phase-shifter 122. A modulating signal X, Y is applied to each of the mixers 120, 121. The output signals of the mixers 120, 121 are applied to a combiner 123 to constitute a modulated signal S. The modulated signal S is either an intermediate frequency signal or a microwave frequency signal. The estimator 124 of the invention includes means 125 for detecting the instantaneous power Pd of the modulated signal S, means 126, 127 for multiplying the detected instantaneous power Pd by each of the modulating signals X, Y, means 128A, 129A for rectifying the signals produced, which rectifier means can be diodes, for example, and means 130, 131 for integrating the rectified signals, the integrator means being followed by subtractors 160, 161 supplying respective signals E1 and E2 proportional to the amplitude of the modulator imbalance. The subtractors 160 and 161 receive respective analog reference signals REF1 and REF2, which are usually identical (REF1=REF2), the output signals of the integrator means 130, 131 being subtracted from REF1 and REF2, respectively. For example, REF1 and REF2 are equal to twice the mean amplitude of the modulated signal S. The subtractors 160 and 161 generally need to be used only for the purposes of correcting imbalance. The signals E1 and E2 are respectively applied to the multipliers 150 and 151 on the input side of the mixers 121 and 120, respectively, to correct the offsets introduced by these mixers.

Although there are already many kinds of compensation circuits and methods, it is still a goal of research to propose newer and better solutions to the I/Q mismatch problem in transmitter systems.

SUMMARY OF THE INVENTION

The object of the present invention is to provide a new method and apparatus of I/Q channel mismatch and local leakage calibration in a transmitter.

The present invention provides a method for I/Q mismatch calibration of a transmitter. First, a discrete-time signal is generated. Next, a corrected signal based on the discrete-time signal and a set of correction parameters A_(p) and B_(p) is obtained. Next, the corrected signal is converted to an analog corrected signal. Next, I/Q modulation is applied to the analog corrected signal and outputs a modulated signal. Next, a first desired component measure and a first image component measure with a first set of the correction parameters A_(p) and B_(p) are obtained from the modulated signal. Next, a second desired component measure and a second image component measure with a second set of the correction parameters A_(p) and B_(p) are obtained from the modulated signal. Next, a third desired component measure and a third image component measure with a third set of the correction parameters A_(p) and B_(p) are obtained from the modulated signal. Next, a fourth and fifth set of correction parameters A_(p) and B_(p) are obtained based on the first, the second, and the third desired component measures as well as the first, the second, and the third image component measures. Next, a fourth desired component measure and a fourth image component measure with the fourth set of correction parameters A_(p) and B_(p) and a fifth desired component measure and a fifth image component measure with the fifth set of correction parameters A_(p) and B_(p) are obtained from the modulated signal. Finally, a final set of the correction parameters A_(p) and B_(p) are selected from the fourth and fifth sets of correction parameters.

In addition, the method for I/Q mismatch calibration of a transmitter according to the present invention further includes the following steps. A DC compensation parameter γ_(p) is added while obtaining the corrected signal x_(c)[n] such that x_(c)[n]=A_(p)·(x[n]+γ_(p))+B_(p)·(x[n]+γ_(p))*. Next, a first local leakage component measure is obtained from the modulated signal with the final set of parameters A_(p) and B_(p), and the first parameter γ_(p)=ζ₁. A second local leakage component measure is then obtained from the modulated signal with the final set of parameters A_(p) and B_(p), and the second parameter γ_(p)=ζ₂. Next, a third local leakage component measure is obtained from the modulated signal with the final set of parameters A_(p) and B_(p), and the third parameter γ_(p)=jζ₁. A fourth local leakage component measure is then obtaining from the modulated signal with the final set of parameters A_(p) and B_(p), and the fourth parameter γ_(p)=jζ₂. Next, a fifth local leakage component measure is obtaining from the modulated signal with the final set of parameters A_(p) and B_(p), and the fifth parameter γ_(p)=0. Finally, a final DC compensation parameter is obtained based on the first local leakage component measure, the second local leakage component measure, the third local leakage component measure, the fourth local leakage component measure and the fifth local leakage component measure.

In addition, the present invention further provides an apparatus for I/Q mismatch calibration of a transmitter. A signal generator generates a discrete-time signal. A correction module receives the discrete-time signal and obtains a corrected signal based on the test signal and a set of correction parameters A_(p) and B_(p). A first and second D/A converter converts the corrected signal to an analog signal, wherein the first D/A converter converts the real part of the corrected signal to the real part of the analog signal, and the second D/A converter converts the imaginary part of the corrected signal to the imaginary part of the analog signal. A modulator applies I/Q modulation to the analog signal, and outputs a modulated signal. A measurer implements the steps of obtaining a first desired component measure and a first image component measure from the modulated signal with a first set of the correction parameters A_(p) and B_(p), obtaining a second desired component measure and a second image component measure from the modulated signal with a second set of the correction parameters A_(p) and B_(p), obtaining a third desired component measure and a third image component measure from the modulated signal with a third set of the correction parameters A_(p) and B_(p), obtaining a fourth desired component measure and a fourth image component measure from the modulated signal with a fourth set of correction parameters A_(p) and B_(p) and obtaining a fifth desired component measure and a fifth image component measure from the modulated signal with a fifth set of correction parameters A_(p) and B_(p). A processor implements the steps of obtaining the fourth and fifth sets of correction parameters A_(p) and B_(p) based on the first, the second, and the third desired component measures as well as the first, the second, and the third image component measures, and selecting a final set of correction parameters A_(p) and B_(p) from the fourth and fifth sets of correction parameters.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will become more fully understood from the detailed description given hereinbelow and the accompanying drawings, given by way of illustration only and thus not intended to be limitative of the present invention.

FIG. 1 shows a quadrature modulator controlled by an imbalance estimator as disclosed in U.S. Pat. No. 6,208,698.

FIG. 2 shows a setup of TX I/Q mismatch and carrier local leakage calibration for 802.11g of a transmitter.

FIGS. 3A-3C are flowcharts of a method for I/Q mismatch and local leakage calibration of a receiver according to one embodiment of the invention.

DETAILED DESCRIPTION OF THE INVENTION

In the I/Q channel transmission architecture, there are two major sources of I/Q mismatch. One is the I/Q modulation and the other is the I/Q transmission band filtering. The mismatch resulting from the I/Q modulation is modeled as a constant over the useful signal bandwidth, which is called systematic I/Q mismatch while the I/Q transmission band filtering causes frequency-dependent I/Q mismatch. The present invention concentrates on the systematic I/Q mismatch. The mathematical theories of the invention will be described in the following.

The baseband signals for the I and Q channels before the I/Q modulation are represented as x_(I)(t) and x_(Q)(t) respectively. Assuming that an amplitude mismatch α, a phase mismatch θ, and zero local leakage results from the I/Q modulation, the equivalent baseband signals y_(I)(t) and y_(Q)(t) for the I and Q channels after the I/Q modulation are, without loss of generality, given by the equations: y _(I)(t)=x _(I)(t)·(1+α)·cos(θ/2)+x _(Q)(t)·(1−α)·sin(θ/2) y _(Q)(t)=x _(I)(t)·(1+α)·sin(θ/2)+x _(Q)(t)·(1−α)·cos(θ/2)  (7)

The combined baseband signal x(t) before the I/Q modulation can be expressed as x_(I)(t)+jx_(Q)(t) while the combined equivalent baseband signal y(t) after the I/Q modulation can be expressed as y_(I)(t)+jy_(Q)(t). Thus, from the equation (7), the I/Q mismatch model can be expressed by the following equation: y(t)=A·x(t)+B·x*(t)  (8) where A and B are complex numbers with

$\begin{matrix} {{A = {\frac{1}{2}\left\{ {{\left( {1 + \alpha} \right){\mathbb{e}}^{j\;{\theta/2}}} + {\left( {1 - \alpha} \right){\mathbb{e}}^{{- j}\;{\theta/2}}}} \right\}}}{B = {\frac{1}{2}\left\{ {{\left( {1 + \alpha} \right){\mathbb{e}}^{j\;{\theta/2}}} - {\left( {1 - \alpha} \right){\mathbb{e}}^{{- j}\;{\theta/2}}}} \right\}}}} & (9) \end{matrix}$ Eq. (8) is for the case of zero local leakage. When the local leakage γ is further considered, the I/Q mismatch model becomes y(t)=A·x(t)+B·x*(t)+γ  (10)

One of solutions to the I/Q mismatch compensation is pre-distortion of the baseband signal x(t) before the I/Q modulation. If the signal is pre-distorted properly, the signal after the I/Q modulation will be free from the I/Q mismatch since the I/Q mismatch and the pre-distortion cancel out. Assuming that the baseband signal x(t) is pre-distorted by a function D, the pre-distorted signal z(t) is expressed as D(x(t)). The ideally calibrated output signal should be:

$\begin{matrix} {{y(t)} = {{A \cdot {D\left( {x(t)} \right)}} + {B \cdot \left( {D\left( {x(t)} \right)} \right)^{*}} + \gamma}} & {{~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~}(11)} \\ {= {C \cdot {x(t)}}} & {(12)} \end{matrix}$ where C is a constant.

Accordingly, the I/Q mismatch compensation can be achieved by identifying the function D to satisfy the equation (12). One of solutions for the function D is a function with parameters A_(p) and B_(p) and is given by: D(x(t))=A _(p) ·x(t)+B _(p) ·x*(t)   (13) As will be shown, Eq. (13) will be implemented via digital signal process. The above is for the case of zero compensation of local leakage. When considering the compensation of the local leakage of the IQ modulation, the function D becomes D(x(t))=A _(p)·(x(t)+γ_(p))+B _(p)·(x(t)+γ_(p))*   (14) By substituting the equation (14) into the equation (11), we derive:

$\begin{matrix} \begin{matrix} {{y(t)} = {{\left( {{A \cdot A_{p}} + {B \cdot B_{p}^{*}}} \right) \cdot {x(t)}} + {\left( {{A \cdot B_{p}} + {B \cdot A_{p}^{*}}} \right) \cdot {x^{*}(t)}} +}} \\ {{\left( {{A \cdot A_{p}} + {B \cdot B_{p}^{*}}} \right) \cdot \gamma_{p}} + {\left( {{A \cdot B_{p}} + {B \cdot A_{p}^{*}}} \right) \cdot \gamma_{p}^{*}} + \gamma} \end{matrix} & (15) \end{matrix}$ The equation (12) is satisfied if

$\begin{matrix} \left\{ \begin{matrix} {\left( {{A \cdot B_{p}} + {B \cdot A_{p}^{*}}} \right) = 0} \\ {\gamma_{p} = {- \gamma}} \end{matrix} \right. & (16) \end{matrix}$

Therefore, the goal of the calibration method in the invention is to identify the coefficients A_(p), B_(p), and γ_(p) in the equation (16).

FIG. 2 shows a setup of TX I/Q channel mismatch and carrier local leakage calibration for 802.11 g of a transmitter according to the embodiment of the present invention. First, a discrete-time signal x[n]=x(nT_(s))=e^(j2πf) ^(T) ^(nT) ^(s) is generated by a signal generator 30, where x(t)=e^(j2πf) ^(T) ^(t) represents a single tone signal, f_(T) is a real number representing the frequency in Hz of the signal x(t) and T_(S) is the sampling period. Next, the signal x[n] is passed to an I/Q correction module 31 including local leakage correction module 32 and gain/phase correction module 33, then the corrected signal x_(c)[n] output by the I/Q correction module 31 is fed to a pair of D/A converters, the first D/A converter 34A and the second D/A converter 34B, which convert the corrected signal x_(c)[n] to an analog signal x_(c)(t). To speak more specifically, the first D/A converter 34A converts the real part of the corrected signal x_(c)[n] to the real part of the analog signal x_(c)(t) and the second D/A converter 34B converts the imaginary part of the corrected signal x_(c)[n] to the imaginary part of the analog signal x_(c)(t). The analog signal x_(c)(t) is applied to an I/Q modulator, which increase the central frequency of the analog signal x_(c)(t) by f_(c) Hz and outputting a modulated signal x_(m)(t), wherein f_(c) is a preset real number. The modulated signal x_(m)(t) is then monitored by a spectrum analyzer 36, or other equipment that can monitor the signal spectrum, to obtain the intensities of the frequency components of the modulated signal x_(m)(t) at f_(c)+f_(T), f_(c)−f_(T), and f_(c) Hz, denoted by W(f_(T)), W(−f_(T)), and L, respectively. The W(f_(T)), W(−f_(T)), and L are actually values indicative of the power of the equivalent baseband signal of x_(m)(t) at f_(T), −f_(T), and 0 Hz, respectively. From the point of view of the I/Q correction, the term W(f_(T)) represents a desired component, which will contain all the energy of the analog corrected signal x_(c)(t) if the I/Q mismatch is ideally compensated. And the term W(−f_(T)) represents the image component, which is the additional interfering component at the image frequency due to the I/Q imbalance. And the term L represents the local leakage component, which is the baseband-equivalent DC component incurred by the carrier leakage. According to the obtained data W(f_(T)), W(−f_(T)), and L, the personal computer (PC) 38 computes the parameters required by the I/Q correction module 31 to minimize the impacts due to the I/Q channel mismatch and local leakage.

FIGS. 3A-3C are flowcharts of a method for I/Q channel mismatch and local leakage calibration of a transmitter according to one embodiment of the invention.

First, a discrete-time signal x[n]=x(nT_(s))=e^(j2πf) ^(T) ^(nT) ^(s) is generated by the signal generator 30 (S1), wherein x(t)=e^(j2πf) ^(T) ^(t) represents a single tone signal, f_(T) is a real number representing the frequency of the signal x(t) and T_(s) is the sampling period. Next, in the present invention, the I/Q gain and phase errors are first calibrated. After the I/Q gain and phase errors are corrected, the local leakage calibration procedure is then performed.

Let x[n]=x_(R)[n]+jx_(I)[n] be the input signal of the I/Q correction module 31, which outputs a corrected signal x_(c)[n] based on the discrete-time signal x[n] and a set of correction parameters (A_(p), B_(p), γ_(p)) (S2). To speak more specifically, the output of the I/Q correction module 31 is expressed by x_(c)[n]=A_(p)·(x[n]+γ_(P))+B_(P)·(x[n]+γ_(P))*

First, the corrected signal x_(c)[n] output by the I/Q correction module 31 is fed to a pair of D/A converters 34, which then convert the corrected signal x_(c)[n] to an analog signal x_(c)(t). To speak more specifically, the first D/A converter converts the real part of the corrected signal x_(c)[n] to the real part of the analog signal x_(c)(t) and the second D/A converter converts the imaginary part of the corrected signal x_(c)[n] to the imaginary part of the analog signal x_(c)(t). The analog signal x_(c)(t) is then applied to an I/Q modulator, which increases the central frequency of the analog signal x_(c)(t) by f_(c) Hz and outputs a modulated signal x_(m)(t), wherein f_(c) is a preset real number. The modulated signal x_(m)(t) is then analyzed by the spectrum analyzer 36. In the present invention, a plurality of sets of parameters A_(p), B_(p) and γ_(p) are predetermined, and then the spectrum analyzer 36 analyzes the modulated signal x_(m)(t) with different sets of the parameters A_(p), B_(p) and γ_(p).

In the I/Q mismatch calibration procedure (S3˜S8) according to the embodiment of the present invention, the parameter γ_(p) is zero during the procedure. A first desired component measure W⁽¹⁾(f_(T)) and a first image component measure W⁽¹⁾(−f_(T)) is obtained with a first set of the correction parameters A_(p) and B_(p) (S3). The first set of the correction parameters A_(p) and B_(p) is (A_(p), B_(p))=(a, 0), wherein a is a real number.

When a=1, we have x_(c)(t)=x(t). The equivalent baseband signal of the modulated signal, denoted by y(t), is given by:

$\begin{matrix} \begin{matrix} {{y(t)} = {{A \cdot {x\lbrack n\rbrack}} + {B \cdot {x^{*}\lbrack n\rbrack}} + \gamma}} \\ {= {{A \cdot {\mathbb{e}}^{j\; 2\;\pi\; f_{T^{t}}}} + {B \cdot {\mathbb{e}}^{{- j}\; 2\;\pi\; f_{T^{t}}}} + \gamma}} \end{matrix} & (17) \end{matrix}$ thus, the intensities of the resulting desired tone component W⁽¹⁾(f_(T)), the image component W⁽¹⁾(−f_(T)) associated with the first set of the correction parameters A_(p) and B_(p) are obtained from the spectrum analyzer 36, see the equation (18).

$\begin{matrix} {\begin{matrix} {{W^{(1)}\left( f_{T} \right)} = {\beta_{1}{A}^{2}}} \\ {= {\beta{{{\cos\frac{\theta}{2}} + {j\;\alpha\;\sin\frac{\theta}{2}}}}^{2}}} \\ {= {\beta\left( {{\cos^{2}\frac{\theta}{2}} + {\alpha^{2}\sin^{2}\frac{\theta}{2}}} \right)}} \end{matrix}\begin{matrix} {{W^{(1)}\left( {- f_{T}} \right)} = {\beta_{1}{B}^{2}}} \\ {= {\beta{{{\alpha\;\cos\frac{\theta}{2}} + {j\;\sin\frac{\theta}{2}}}}^{2}}} \\ {= {\beta\left( {{\alpha^{2}\cos^{2}\frac{\theta}{2}} + {\sin^{2}\frac{\theta}{2}}} \right)}} \end{matrix}} & (18) \end{matrix}$ where β₁ is an unknown constant related to the spectrum analyzer 36. Here, the measures W⁽¹⁾(f_(T)) and W⁽¹⁾(−f_(T)) are values indicative of the power of the signal y(t) at f_(T) and −f_(T) Hz, respectively.

In this case, the image-rejection ratio ρ⁽¹⁾ is

$\begin{matrix} {\rho^{(1)} = {\frac{W^{(1)}\left( {- f_{T}} \right)}{W^{(1)}\left( f_{T} \right)} = \frac{{\alpha^{2}\cos^{2}\frac{\theta}{2}} + {\sin^{2}\frac{\theta}{2}}}{{\cos^{2}\frac{\theta}{2}} + {\alpha^{2}\sin^{2}\frac{\theta}{2}}}}} & (19) \end{matrix}$

Next, a second desired component measure W⁽²⁾(f_(T)) and a second image component measure W⁽²⁾(−f_(T)) are obtained from the spectrum analyzer 36 with a second set of the correction parameters A_(p) and B_(p) (S4). The second set of the correction parameters A_(p) and B_(p) is (A_(p), B_(p))=(b, b), wherein ‘b’ is a real number.

${{When}\mspace{14mu} b} = {\frac{1}{2}\text{:}}$ $\begin{matrix} \begin{matrix} {{y(t)} = {{A \cdot \left( {{0.5 \cdot {x(t)}} + {0.5 \cdot {x^{*}(t)}}} \right)} + {B \cdot \left( {{0.5 \cdot {x(t)}} + {0.5 \cdot {x^{*}(t)}}} \right)^{*}} + \gamma}} \\ {= {{\left( {{0.5A} + {0.5B}} \right) \cdot {x(t)}} + {\left( {{0.5A} + {0.5B}} \right) \cdot {x^{*}(t)}} + \gamma}} \\ {= {{\left( {{0.5A} + {0.5B}} \right) \cdot {\mathbb{e}}^{j\; 2\;\pi\; f_{T^{t}}}} + {\left( {{0.5A} + {0.5B}} \right) \cdot {\mathbb{e}}^{{- j}\; 2\;\pi\; f_{T^{t}}}} + \gamma}} \end{matrix} & (20) \end{matrix}$

Thus, the intensities of the resulting desired tone component W⁽²⁾(f_(T)), the image component W⁽²⁾(−f_(T)) are obtained from the spectrum analyzer 36 with the second set of the correction parameters A_(p) and B_(p), see the equation (21).

$\begin{matrix} {\begin{matrix} {{W^{(2)}\left( f_{T} \right)} = {\beta_{2}{{{0.5A} + {0.5B}}}^{2}}} \\ {= {\frac{\beta_{2}}{4}{{{\cos\frac{\theta}{2}} + {j\;\alpha\;\sin\frac{\theta}{2}} + {\alpha\;\cos\frac{\theta}{2}} + {j\;\sin\frac{\theta}{2}}}}^{2}}} \\ {= {\frac{\beta_{2}}{4}{{{\left( {1 + \alpha} \right)\cos\frac{\theta}{2}} + {{j\left( {1 + \alpha} \right)}\sin\frac{\theta}{2}}}}^{2}}} \\ {= {\frac{\beta_{2}}{4}\left( {{\left( {1 + \alpha} \right)^{2}\cos^{2}\frac{\theta}{2}} + {\left( {1 + \alpha} \right)^{2}\sin^{2}\frac{\theta}{2}}} \right)}} \\ {= \frac{{\beta_{2}\left( {1 + \alpha} \right)}^{2}}{4}} \end{matrix}{{W^{(2)}\left( {- f_{T}} \right)} = {{\beta_{2}{{{0.5A} + {0.5B}}}^{2}} = \frac{{\beta_{2}\left( {1 + \alpha} \right)}^{2}}{4}}}} & (21) \end{matrix}$

Here, from the output of the spectrum analyzer 36, a value N=W⁽²⁾(−f_(T))=W⁽²⁾(f_(T)) is proportional to (1+α)².

Next, a third desired component measure W⁽³⁾(f_(T)) and a third image component measure W⁽³⁾(−f_(T)) are obtained from the spectrum analyzer 36 with a third set of the correction parameters A_(p) and B_(p) (S5). The third set of the correction parameters A_(p) and B_(p) is (A_(p), B_(p))=(b, −b).

${{When}\mspace{14mu} b} = {\frac{1}{2}\text{:}}$ $\begin{matrix} \begin{matrix} {{y(t)} = {{A\left( {{0.5 \cdot {x(t)}} - {0.5 \cdot {x^{*}(t)}}} \right)} + {B\left( {{0.5 \cdot {x(t)}} - {0.5 \cdot {x^{*}(t)}}} \right)}^{*} + \gamma}} \\ {= {{\left( {{0.5A} - {0.5B}} \right) \cdot {x(t)}} + {\left( {{{- 0.5}A} + {0.5B}} \right) \cdot {x^{*}(t)}} + \gamma}} \\ {= {{\left( {{0.5A} - {0.5B}} \right) \cdot {\mathbb{e}}^{j\; 2\;\pi\; f_{T^{t}}}} + {\left( {{{- 0.5}A} + {0.5B}} \right) \cdot {\mathbb{e}}^{{- j}\; 2\;\pi\; f_{T^{t}}}} + \gamma}} \end{matrix} & (22) \end{matrix}$

Thus, the intensities of the resulting desired tone component W⁽³⁾(f_(T)), the image component W⁽³⁾(−f_(T)) are obtained from the spectrum analyzer 36 with the third set of the correction parameters A_(p) and B_(p), see the equation (23).

$\begin{matrix} {\begin{matrix} {{W^{(3)}\left( f_{T} \right)} = {\beta_{2}{{{0.5A} - {0.5B}}}^{2}}} \\ {= {\frac{\beta_{2}}{4}{{{\cos\frac{\theta}{2}} + {j\;\alpha\;\sin\frac{\theta}{2}} - {\alpha\;\cos\frac{\theta}{2}} - {j\;\sin\frac{\theta}{2}}}}^{2}}} \\ {= {\frac{\beta_{2}}{4}{{{\left( {1 - \alpha} \right)\cos\frac{\theta}{2}} + {{j\left( {\alpha - 1} \right)}\sin\frac{\theta}{2}}}}^{2}}} \\ {= {\frac{\beta_{2}}{4}\left( {{\left( {1 - \alpha} \right)^{2}\cos^{2}\frac{\theta}{2}} + {\left( {1 - \alpha} \right)^{2}\sin^{2}\frac{\theta}{2}}} \right)}} \\ {= \frac{{\beta_{2}\left( {1 - \alpha} \right)}^{2}}{4}} \end{matrix}{{W^{(3)}\left( {- f_{T}} \right)} = {{\beta_{2}{{{{- 0.5}A} + {0.5B}}}^{2}} = \frac{{\beta_{2}\left( {1 - \alpha} \right)}^{2}}{4}}}} & (23) \end{matrix}$ where β₂ is an unknown constant related to the spectrum analyzer 36. Here, from the output of the spectrum analyzer 36, a value O=W⁽³⁾(−f_(T))=W⁽³⁾(f_(T)) is proportional to (1−α)². Thus, the gain mismatch term α can be estimated by

$\begin{matrix} {{\alpha \approx \hat{\alpha}} = \frac{\sqrt{\frac{N}{O}} - 1}{\sqrt{\frac{N}{O}} + 1}} & (24) \end{matrix}$

Term N can be set to be either N=W⁽²⁾(f_(T)) or N=W⁽²⁾(−f_(T)). Due to the existence of the measurement error, the term N can also set to be N=(W⁽²⁾(f_(T))+W⁽²⁾(−f_(T)))/2 or N=√{square root over (W⁽²⁾(f_(T))·W⁽²⁾(−f_(T)))}{square root over (W⁽²⁾(f_(T))·W⁽²⁾(−f_(T)))}{square root over (W⁽²⁾(f_(T))·W⁽²⁾(−f_(T)))}{square root over (W⁽²⁾(f_(T))·W⁽²⁾(−f_(T)))} to increase the accuracy. Similarly, Term O can be set to be either O=W⁽³⁾(f_(T)) or O=W⁽³⁾(−f_(T)). Due to the existence of the measurement error, the term O can also set to be O=(W⁽³⁾(f_(T))+W⁽³⁾(−f_(T)))/2 or O=√{square root over (W⁽³⁾(f_(T))·W⁽³⁾(−f_(T)))}{square root over (W⁽³⁾(f_(T))·W⁽³⁾(−f_(T)))}{square root over (W⁽³⁾(f_(T))·W⁽³⁾(−f_(T)))}{square root over (W⁽³⁾(f_(T))·W⁽³⁾(−f_(T)))} to increase the accuracy. Substitute α into the equation (19), the values Q and P are obtained, see the equation (25).

$\begin{matrix} {Q = {\frac{{\hat{\alpha}}^{2} - \rho^{(1)}}{\left( {1 + \rho^{(1)}} \right)\left( {{\hat{\alpha}}^{2} - 1} \right)} = {\sin^{2}\frac{\theta}{2}}}} & (25) \\ {P = {{1 - Q} = {{1 - {\sin^{2}\frac{\theta}{2}}} = {\cos^{2}\frac{\theta}{2}}}}} & (26) \end{matrix}$

Next, a fourth and fifth set of correction parameters A_(p) and B_(p) are obtained based on the first desired component measure W⁽¹⁾(f_(T)), the second desired component measure W⁽²⁾(f_(T)), and the third desired component measure W⁽³⁾(f_(T)), and the first image component measure W⁽¹⁾(−f_(T)), the second image component measure W^((2)(−f) _(T)), and the third image component measure W⁽³⁾(−f_(T)) (S6).

The fourth and fifth sets of correction parameters A_(p) and B_(p) are expressed by parameters {circumflex over (α)}, P, and Q. The fourth set of the correction parameters A_(p) and B_(p) is A _(p) =√{square root over (P)}−j{circumflex over (α)}√{square root over (Q)} B _(p) =−{circumflex over (α)}√{square root over (P)}−j√{square root over (Q)}, and the fifth set of the correction parameters A_(p) and B_(p) is A _(p) =√{square root over (P)}+j{circumflex over (α)}√{square root over (Q)} B _(p) =−{circumflex over (α)}√{square root over (P)}+j√{square root over (Q)},

Next, a fourth desired component measure W⁽⁴⁾(f_(T)) and a fourth image component measure W⁽⁴⁾(−f_(T)) are obtained with the fourth set of correction parameters A_(p) and B_(p), and a fifth desired component measure W⁽⁵⁾(f_(T)) and a fifth image component measure W⁽⁵⁾(−f_(T)) are obtained with the fifth set of correction parameters A_(p) and B_(p) (S7).

Last, a final set of the correction parameters A_(p) and B_(p) is determined according to the fourth and fifth sets of correction parameters (S8). In addition, since there is no information about the sign of the phase mismatch θ, a hypothesis testing is performed to determine whether

${{\sin\frac{\theta}{2}} = {{\sqrt{Q}\mspace{14mu}{or}\mspace{14mu}\sin\frac{\theta}{2}} = {- \sqrt{Q}}}},$ however, in either case,

${\cos\frac{\theta}{2}} = \sqrt{1 - Q}$ since

${- \frac{\pi}{2}} < \theta < \frac{\pi}{2}$ in typical conditions.

If the image rejection ratio corresponding to the fourth set of (A_(p), B_(p)) is less than or equal to the image rejection ratio corresponding to the fifth set of (A_(p), B_(p)), the fourth set of (A_(p), B_(p)) is chosen as the final set of (A_(p), B_(p)), otherwise, the fifth set of (A_(p), B_(p)) is chosen as the final set of (A_(p), B_(p)). In addition, the final set of (A_(p), B_(p)) can be further multiplied by a normalization factor

$\chi = {\frac{1}{{AA}_{p} + {BB}_{p}^{*}} = \frac{1}{\left( {1 - {\hat{\alpha}}^{2}} \right) \cdot \left( {1 - {2Q}} \right)}}$ to ensure the unity gain of the correction process.

The steps S6 and S8 are performed by the personal computer 38 shown in FIG. 2 of the present invention.

In addition, it is not the only way to determine the final set of the correction parameters A_(p) and B_(p) through the value of

$\frac{W^{(4)}\left( {- f_{T}} \right)}{W^{(4)}\left( f_{T} \right)}\mspace{14mu}{and}\mspace{14mu}{\frac{W^{(5)}\left( {- f_{T}} \right)}{W^{(5)}\left( f_{T} \right)}.}$ Under the final parameter, the image component would be very small (theoretically it is zero), and therefore would be much less than the image component under the other sets of parameter. From another angle, a small image component will imply a strong desired component, since their sum should be the same in split regardless of the set of parameters. Therefore, under the final parameter, the desired component would be stronger than the desired component under the other set of parameters.

Next, the local leakage calibration is active after the elimination of I/Q channel mismatch errors in step S9˜S13.

In the local leakage calibration procedure according to the embodiment of the present invention, a first local leakage component measure L₁ of the corrected signal is obtained from the spectrum analyzer 36 with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=ζ₁ (S9), where ζ₁ is a real number.

Thus, the carrier leakage term at the spectrum analyzer 36 is L ₁=β₃·(γ_(R) ²+γ_(I) ²+ζ₁ ²+2γ_(R)·ζ₁).  (27) where β₃ is a scaling factor.

Next, a second local leakage component measure L₂ of the corrected signal is obtained from the spectrum analyzer 36 with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=ζ₂ (S10), where ζ₂ is a real number.

Thus, the carrier leakage term at the spectrum analyzer 36 is L ₂=β₃·(γ_(R) ²+γ_(I) ²+ζ₂ ²+2γ_(R)·ζ₂).  (28)

Next, a third local leakage component measure L₃ of the corrected signal is obtained from the spectrum analyzer 36 with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=jζ₁ (S11), where ζ₁ is a real number.

Thus, the carrier leakage term at the spectrum analyzer 36 is L ₃=β₃·(γ_(R) ²+γ_(I) ²+ζ₁ ²+2γ_(I)·ζ₁).  (29)

Next, a fourth local leakage component measure L₄ of the corrected signal is obtained from the spectrum analyzer 36 with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=jζ₂ (S12), where ζ₂ is another real number.

Thus, the carrier leakage term at the spectrum analyzer 36 is L ₄=β₃·(γ_(R) ²+γ_(I) ²+ζ₂ ²+2γ_(I)·ζ₂).  (30)

Next, a fifth local leakage component measure L₅ of the corrected signal is obtained from the spectrum analyzer 36 with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=0 (S13).

Thus, the carrier leakage term at the spectrum analyzer 36 is L ₅=β₃·(γ_(R) ²+γ_(I) ²).  (31)

Thus, the parameters γ_(R) and γ_(I) are obtained based on the first local leakage component measure L₁, the second local leakage component measure L₂, the third local leakage component measure L₃, the fourth local leakage component measure L₄ and the fifth local leakage component measure L₅ (S14), see the equations (32) and (33).

$\begin{matrix} {{\hat{\gamma}}_{R} = {\frac{1}{2}\mspace{11mu}\frac{{\zeta_{2}^{2}\left( {L_{1} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{2} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{2} - L_{5}} \right)} - {\zeta_{2}\left( {L_{1} - L_{5}} \right)}}}} & (32) \\ {{\hat{\gamma}}_{I} = {{- \frac{1}{2}}\mspace{11mu}\frac{{\zeta_{2}^{2}\left( {L_{3} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{4} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{4} - L_{5}} \right)} - {\zeta_{2}\left( {L_{3} - L_{5}} \right)}}}} & (33) \end{matrix}$

Finally, according to the step S14, the final DC compensation parameter γ_(p,final) is obtained by

$\begin{matrix} \begin{matrix} {\gamma_{p,{final}} = {{{- \frac{1}{2}} \cdot \frac{{\zeta_{2}^{2}\left( {L_{1} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{2} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{2} - L_{5}} \right)} - {\zeta_{2}\left( {L_{1} - L_{5}} \right)}}} -}} \\ {j{\frac{1}{2} \cdot \frac{{\zeta_{2}^{2}\left( {L_{3} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{4} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{4} - L_{5}} \right)} - {\zeta_{2}\left( {L_{3} - L_{5}} \right)}}}} \end{matrix} & (34) \end{matrix}$

In addition, the step S14 is performed by the personal computer 38 shown in FIG. 2 of the present invention.

The previously described method can be applied to a transceiver module, e.g., an IEEE 802.11 compliant wireless LAN transceiver module, involving I/Q demodulation. Wireless LANs based on the IEEE 802.11 standard have achieved wide customer acceptance in the enterprise environment. They are expected to continue to expand in popularity and become ubiquitous communication systems even in private and public places.

Thus, after the calibration procedure, the I/Q correction module compensates the I/Q channel mismatch in the signal from the baseband transmitter by the finally determined coefficients A_(p), B_(p) and γ_(p). The previously described calibration procedure may be implemented upon production stage of the transceiver.

In conclusion, the present invention provides a new method and apparatus for transmitter I/Q mismatch calibration, especially suitable for an IEEE 802.11g compliant WLAN transceiver module. The compensation of the I/Q mismatch is achieved by a signal pre-distortion prior to the I/Q modulation. The coefficients of the pre-distortion function are determined by a calibration procedure. After performing the procedure disclosed by the embodiment of the present invention, the coefficients are accurate enough for the transmitted signal to be free from the I/Q mismatch.

In addition, the present invention can be applied to any kind of transmitter having an I/Q modulator, and the OFDM transmitter is just one example, thus it is not limited to the scope of the present invention.

The foregoing description of the preferred embodiments of this invention has been presented for purposes of illustration and description. Obvious modifications or variations are possible in light of the above teaching. The embodiments were chosen and described to provide the best illustration of the principles of this invention and its practical application to thereby enable those skilled in the art to utilize the invention in various embodiments and with various modifications as are suited to the particular use contemplated. All such modifications and variations are within the scope of the present invention as determined by the appended claims when interpreted in accordance with the breadth to which they are fairly, legally, and equitably entitled. 

1. A method for in-phase and quadrature mismatch calibration of a transmitter, comprising the following steps: generating a discrete-time signal x[n]=x(n·T_(s)), wherein x(t)=e^(j2πf) ^(T) ^(t) and f_(T) and T_(s) are real numbers; obtaining a corrected signal x_(c)[n] based on the signal x[n] and a set of correction parameters A_(p) and B_(p), wherein x _(c) [n]=A _(p) ·x[n]+B _(p) ·x*[n]; converting the corrected signal x_(c)[n] to an analog corrected signal x_(c)(t); applying in-phase and quadrature modulation to the analog corrected signal x_(c)(t) and outputting a modulated signal x_(m)(t); obtaining a first desired component measure W⁽¹⁾(f_(T)) and a first image component measure W⁽¹⁾(−f_(T)) from the modulated signal x_(m)(t) with a first set of the correction parameters A_(p) and B_(p); obtaining a second desired component measure W⁽²⁾(f_(T)) and a second image component measure W⁽²⁾(−f_(T)) from the modulated signal x_(m)(t) with a second set of the correction parameters A_(p) and B_(p); obtaining a third desired component measure W⁽³⁾(f_(T)) and a third image component measure W⁽³⁾(−f_(T)) from the modulated signal x_(m)(t) with a third set of the correction parameters A_(p) and B_(p); obtaining a fourth and fifth set of correction parameters A_(p) and B_(p) based on the first, the second, and the third desired component measures as well as the first, the second, and the third image component measures; obtaining a fourth desired component measure W⁽⁴⁾(f_(T)) and a fourth image component measure W⁽⁴⁾(−f_(T)) from the modulated signal x_(m)(t) with the fourth set of correction parameters A_(p) and B_(p); obtaining a fifth desired component measure W⁽⁵⁾(f_(T)) and a fifth image component measure W⁽⁵⁾(−f_(T)) from the modulated signal x_(m)(t) with the fifth set of correction parameters A_(p) and B_(p); and obtaining a final set of the correction parameters A_(p) and B_(p) from the fourth and fifth sets of correction parameters.
 2. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 1, wherein the first set of correction parameters (A_(p), B_(p))=(a,0), the second set of correction parameters (A_(p),B_(p))=(b,b), and the third set of correction parameters (A_(p),B_(p))=(b,−b), where a and b are real numbers.
 3. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 2, wherein the parameter a is 1 and the parameter b is ½.
 4. The method for mismatch calibration of a transmitter as claimed in claim 1, wherein the fourth set of correction parameters (A_(p),B_(p)) are obtained by A _(p) =√{square root over (P)}−j{circumflex over (α)}√{square root over (Q)} B _(p) =−{circumflex over (α)}√{square root over (P)}−j√{square root over (Q)} and the fifth set of correction parameters (A_(p),B_(p)) are obtained by A _(p) =√{square root over (P)}+j{circumflex over (α)}√{square root over (Q)} B _(p) =−{circumflex over (α)}√{square root over (P)}+j√{square root over (Q)} where ${{\alpha \approx \hat{\alpha}} = \frac{\sqrt{\frac{N}{O}} - 1}{\sqrt{\frac{N}{O}} + 1}},{N = {\left( {{W^{(2)}\left( f_{T} \right)} + {W^{(2)}\left( {- f_{T}} \right)}} \right)/2}},{O = {\left( {{W^{(3)}\left( f_{T} \right)} + {W^{(3)}\left( {- f_{T}} \right)}} \right)/2}},{Q = \frac{{\hat{\alpha}}^{2} - \rho^{(1)}}{\left( {1 + \rho^{(1)}} \right)\left( {{\hat{\alpha}}^{2} - 1} \right)}},{P = {1 - Q}},{\rho^{(1)} = {\frac{W^{(1)}\left( {- f_{T}} \right)}{W^{(1)}\left( f_{T} \right)}.}}$
 5. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 1, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a function of W⁽⁴⁾(−f_(T)) is less than the function of W⁽⁵⁾(−f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 6. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 5, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a value of W⁽⁴⁾(−f_(T)) is less than a value of W⁽⁵⁾(−f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 7. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 1, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a function of W⁽⁴⁾(f_(T)) is greater than the function of W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 8. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 7, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a value of W⁽⁴⁾(f_(T)) is greater than W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 9. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 1, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a function of W⁽⁴⁾(−f_(T)) and W⁽⁴⁾(f_(T)) is less than the function of W⁽⁵⁾(−f_(T)) and W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 10. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 9, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if W⁽⁴⁾(−f_(T))/W⁽⁴⁾(f_(T)) is less than W⁽⁵⁾(−f_(T))/W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 11. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 1, further comprising the following steps: further adding an DC compensation parameter γ_(p) while obtaining the corrected signal x_(c)[n] such that x _(c) [n]=A _(p)·(x[n]+γ _(p))+B _(p)·(x[n]+γ _(p))*; obtaining a first local leakage component measure L₁ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=ζ₁, where ζ₁ is a real number; obtaining a second local leakage component measure L₂ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=ζ₂, where ζ₂ is a real number; obtaining a third local leakage component measure L₃ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=jζ₁; obtaining a fourth local leakage component measure L₄ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=jζ₂; obtaining a fifth local leakage component measure L₅ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=0; and obtaining a final DC compensation parameter γ_(p,final) based on the first local leakage component measure L₁, the second local leakage component measure L₂, the third local leakage component measure L₃, the fourth local leakage component measure L₄ and the fifth local leakage component measure L₅.
 12. The method for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 11, wherein the final DC compensation parameter γ_(p,final) is obtained by $\begin{matrix} {\gamma_{p,{final}} = {{{- \frac{1}{2}} \cdot \frac{{\zeta_{2}^{2}\left( {L_{1} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{2} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{2} - L_{5}} \right)} - {\zeta_{2}\left( {L_{1} - L_{5}} \right)}}} -}} \\ {j{\frac{1}{2} \cdot \frac{{\zeta_{2}^{2}\left( {L_{3} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{4} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{4} - L_{5}} \right)} - {\zeta_{2}\left( {L_{3} - L_{5}} \right)}}}} \end{matrix}.$
 13. An apparatus for in-phase and quadrature mismatch calibration of a transmitter, comprising: a signal generator for generating a discrete-time signal x[n]=x(n·T_(s)), wherein x(t)=e^(j2πf) ^(T) ^(t) and f_(▭) and T_(s) are real numbers; a correction module for receiving the discrete-time signal x[n] and obtaining a corrected signal x_(c)[n] based on the test signal x[n] and a set of correction parameters A_(p) and B_(p), wherein x _(c) [n]=A _(p) ·x[n]+B _(p) ·x*[n]; a first and second D/A converter converting the corrected signal x_(c)[n] to an analog signal x_(c)(t), wherein the first D/A converter converts the real part of the corrected signal to a real part of the analog signal, and the second D/A converter converts the imaginary part of the corrected signal to an imaginary part of the analog signal; a modulator applying in-phase and quadrature modulation to the analog signal x_(c)(t), and outputting a modulated signal x_(m)(t); a measurer configured to: obtain a first desired component measure W⁽¹⁾(f_(T)) and a first image component measure W⁽¹⁾(−f_(T)) from the modulated signal x_(m)(t) with a first set of the correction parameters A_(p) and B_(p); obtain a second desired component measure W⁽²⁾(f_(T)) and a second image component measure W⁽²⁾(−f_(T)) from the modulated signal x_(m)(t) with a second set of the correction parameters A_(p) and B_(p); obtain a third desired component measure W⁽³⁾(f_(T)) and a third image component measure W⁽³⁾(−f_(T)) from the modulated signal x_(m)(t) with a third set of the correction parameters A_(p) and B_(p); obtain a fourth desired component measure W⁽⁴⁾(f_(T)) and a fourth image component measure W⁽⁴⁾(−f_(T)) from the modulated signal x_(m)(t) with a fourth set of correction parameters A_(p) and B_(p); and obtain a fifth desired component measure W⁽⁵⁾(f_(T)) and a fifth image component measure W⁽⁵⁾(−f_(T)) from the modulated signal x_(m)(t) with a fifth set of correction parameters A_(p) and B_(p); and a processor configured to: obtain the fourth and fifth sets of correction parameters A_(p) and B_(p) based on the first, the second, and the third desired component measures as well as the first, the second, and the third image component measures; and choose a final set of correction parameters A_(p) and B_(p) from the fourth and fifth sets of correction parameters.
 14. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 13, wherein the first set of correction parameters (A_(p),B_(p))=(a,0), the second set of correction parameters (A_(p),B_(p))=(b,b), and the third set of correction parameters (A_(p),B_(p))=(b,−b), where a and b are real numbers.
 15. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 14, wherein the parameter a is 1 and the parameter b is ½.
 16. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 13, wherein the fourth set of correction parameters (A_(p),B_(p)) are obtained by A _(p) =√{square root over (P)}−j{circumflex over (α)}√{square root over (Q)} B _(p)=−{circumflex over (α)}√{square root over (P)}−j√{square root over (Q)} and the fifth set of correction parameters (A_(p),B_(p)) are obtained by A _(p) =√{square root over (P)}+j{circumflex over (α)}√{square root over (Q)} B _(p) =−{circumflex over (α)}√{square root over (P)}+j√{square root over (Q)} where ${{\alpha \approx \hat{\alpha}} = \frac{\sqrt{N/O} - 1}{\sqrt{N/O} + 1}},{N = {\left( {{W^{(2)}\left( f_{T} \right)} + {W^{(2)}\left( {- f_{T}} \right)}} \right)/2}},{O = {\left( {{W^{(3)}\left( f_{T} \right)} + {W^{(3)}\left( {- f_{T}} \right)}} \right)/2}},{Q = \frac{{\hat{\alpha}}^{2} - \rho^{(1)}}{\left( {1 + \rho^{(1)}} \right)\left( {{\hat{\alpha}}^{2} - 1} \right)}},{P = {1 - Q}},{\rho^{(1)} = {\frac{W^{(1)}\left( {- f_{T}} \right)}{W^{(1)}\left( f_{T} \right)}.}}$
 17. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 13, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a function of W⁽⁴⁾(−f_(T)) is less than the function of W⁽⁵⁾(−f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 18. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 17, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a value of W⁽⁴⁾(−f_(T)) is less than a value of W⁽⁵⁾(−f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 19. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 13, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a function of W⁽⁴⁾(−f_(T)) is greater than the function of W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 20. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 19, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a value of W⁽⁴⁾(f_(T)) is greater than a value of W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 21. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 13, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if a function of W⁽⁴⁾(−f_(T)) and W⁽⁴⁾(f_(T)) is less than the function of W⁽⁵⁾(−f_(T)) and W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 22. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 21, wherein the final set of correction parameters (A_(p),B_(p)) is set to be the fourth set of correction parameters if W⁽⁴⁾(−f_(T))/W⁽⁴⁾(f_(T)) is less than W⁽⁵⁾(−f_(T))/W⁽⁵⁾(f_(T)), otherwise the final set of correction parameters (A_(p),B_(p)) is set to be the fifth set of correction parameters.
 23. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 13, wherein the processor further configured to: further add an DC compensation parameter γ_(p) while obtaining the corrected signal x_(c)[n] such that x _(c) [n]=A _(p)·(x[n]+γ _(p))+B _(p)·(x[n]+γ _(p))*; obtain a first local leakage component measure L₁ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=ζ₁, where ζ₁ is a real number; obtain a second local leakage component measure L₂ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)ζ₂, where ζ₂ is a real number; obtain a third local leakage component measure L₃ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=jζ₁; obtain a fourth local leakage component measure L₄ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=jζ₂; obtain a fifth local leakage component measure L₅ from the modulated signal x_(m)(t) with the final set of parameters A_(p) and B_(p), and the parameter γ_(p)=0; and obtain a final DC compensation parameter γ_(p,final) based on the first local leakage component measure L₁, the second local leakage component measure L₂, the third local leakage component measure L₃, the fourth local leakage component measure L₄ and the fifth local leakage component measure L₅.
 24. The apparatus for in-phase and quadrature mismatch calibration of a transmitter as claimed in claim 23, wherein the final DC compensation parameter γ_(p,final) is obtained by $\gamma_{p,{final}} = {{{- \frac{1}{2}} \cdot \frac{{\zeta_{2}^{2}\left( {L_{1} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{2} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{2} - L_{5}} \right)} - {\zeta_{2}\left( {L_{1} - L_{5}} \right)}}} - {j{\frac{1}{2} \cdot {\frac{{\zeta_{2}^{2}\left( {L_{3} - L_{5}} \right)} - {\zeta_{1}^{2}\left( {L_{4} - L_{5}} \right)}}{{\zeta_{1}\left( {L_{4} - L_{5}} \right)} - {\zeta_{2}\left( {L_{3} - L_{5}} \right)}}.}}}}$ 